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In this paper we present a novel space radiation shielding approach using various 
material lay-ups, called “Graded-Z” shielding, which could optimize cost, weight, and safety 
while mitigating the radiation exposures from the trapped radiation and solar proton 
environments, as well as the galactic cosmic radiation (GCR) environment, to humans and 
electronics. In addition, a validation and verification (V&V) was performed using two 
different high energy particle transport/dose codes (MCNPX & HZETRN). Inherently, we 
know that materials having high-hydrogen content are very good space radiation shielding 
materials. Graded-Z material lay-ups are very good trapped electron mitigators for medium 
earth orbit (MEO) and geostationary earth orbit (GEO). In addition, secondary particles, 
namely neutrons, are produced as the primary particles penetrate a spacecraft, which can 
have deleterious effects to both humans and electronics. The use of “dopants,” such as 
beryllium, boron, and lithium, impregnated in other shielding materials provides a means of 
absorbing the secondary neutrons. Several examples of optimized Graded-Z shielding lay- 
ups that include the use of composite materials are presented and discussed in detail. This 
parametric shielding study is an extension of some earlier pioneering work we (William 
Atwell and Kristina Rojdev) performed in 2004 1 and 2009 2 . 

Nomenclature 

GCR = galactic cosmic radiation 

GEO = geostationary earth orbit 

HZETRN = High Z and Energy Transport code 

MCNPX = Monte Carlo Nuclear Particle transport computer code 

MEO = medium earth orbit 

V & V = validation & verification 


I. Introduction 

G RADED-Z shielding is a laminate of several materials with different Z values (atomic numbers) designed to 
protect against ionizing radiation. Compared to single-material shielding, the same mass of Graded-Z shielding 
has been shown to reduce electron penetration over 60%. 3 Graded-Z shielding is commonly used in satellite- 
based particle detectors and offers several benefits: protection from radiation damage, reduction of background noise 
for detectors, and lower mass compared to single-material shielding. Designs vary, but typically involve a gradient 
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from high-Z (usually tantalum [Ta]) through successively lower-Z elements such as tin, steel, and copper, usually 
ending with aluminum. Sometimes even lighter materials such as polypropylene and boron carbide are used. 4,5 In a 
typical Graded-Z shield, the high-Z layer effectively scatters protons and electrons. It also absorbs gamma rays, 
which produces x-ray fluorescence. Each subsequent layer absorbs the X-ray fluorescence of the previous material, 
eventually reducing the energy to a suitable level. Each decrease in energy produces bremsstrahlung and Auger 
electrons, which are below the detector's energy threshold. Some designs also include an outer layer of aluminum, 
which may simply be the outer structure of a spacecraft or satellite. Fan, et al. 6 reports “Shielding for space 
microelectronics needs to provide an acceptable dose rate with minimum shield mass. The analysis presented here 
shows that the best approach is, in general, to use a graded-Z shield, with a high-Z layer sandwiched between two 
low-Z materials. A Graded-Z shield is shown to reduce the electron dose rate by more than sixty percent over a 
single-material shield of the same areal density. For protons, the optimal shield would consist of a single, low-Z 
material layer. However, it is shown that a Graded-Z shield is nearly as effective as a single-material shield, as long 
as a low-Z layer is located adjacent to the microelectronics. A specific shield design depends upon the details of the 
radiation environment, system model, design margins/levels, compatibility of shield materials, etc.” 

In this paper we investigate several novel Graded-Z materials that are compared with baseline materials: 
aluminum (Al), high density polyethylene (HDPE), and water (H 2 0). In addition, we compare the results using two 
different computational computer codes, NASA HZETRN 7 and MCNPX 8 to validate and verify our results. 


II. Study Assumptions 


A. Graded-Z configurations 

The shielding configurations considered in this paper were two-layer shields with typical shielding thicknesses 
and very large thickness configurations. The two layer shielding configurations were Al-HDPE, Ta-HDPE, and W- 
HDPE. The thicknesses of Al, Ta, and W were 150 mils, which correspond to 1.0287, 6.3237, and 7.3343 g/cm 2 , 
respectively, and the HDPE thicknesses were 5 and 10 g/cm 2 . The very large thickness configurations consisted of 
single layer H 2 0 and two-layer H 2 0-SWX. SWX, specifically SWX-210, is a boron-doped hydrocarbon having a 
density of 1.19 g/cm 3 . The water thicknesses ranged from 100-300 g/cm 2 , and the SWX thicknesses were 10 and 20 
g/cm 2 . 


B. Environments 

For deep-space exploratory missions, the radiation sources are due to the naturally-occurring galactic cosmic 
radiation (GCR) and from solar particle events (SPEs). The radiation exposure analyses discussed later is based on 
the GCR environment during solar minimum (epoch 1977), solar maximum (epoch 1982), and the classic series of 
solar proton events that occurred 19-24 October 1989. The October 1989 SPEs were four events that were combined 
using the Band fit. 9,10 For all three environment cases, the energy spectrum was restricted to 2 GeV due to 
limitations in the MCNPX software, as shown in the differential energy spectrum below (Figure 1). 
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Figure 1. The Band fit differential energy spectrum for the 19-24 October 1989 SPEs. 



C. High Energy Particle Transport/Dose Codes 

Two high energy particle transport/dose codes were used in the computations, HZETRN 7 and MCNPX 8 , and the 
results were compared as a check on the verification and validation of the absorbed doses. 

The HZETRN code was develop at NASA Langley Research Center in 1995. The high-charge-and energy (HZE) 
transport computer program HZETRN was developed to address the problems of free-space radiation transport and 
shielding. The HZETRN program is intended specifically for the design engineer who is interested in obtaining fast 
and accurate dosimetric information for the design and construction of space modules and devices. The program is 
based on a one-dimensional space-marching formulation of the Boltzmann transport equation with a straight-ahead 
approximation. The effect of the long-range Coulomb force and electron interaction is treated as a continuous 
slo wing-down process. Atomic (electronic) stopping power coefficients with energies above a few MeV are 
calculated by using Bethe' s theory including Bragg' s rule, Ziegler' s shell corrections, and effective charge. Nuclear 
absorption cross sections are obtained from fits to quantum calculations and total cross sections are obtained with a 
Ramsauer formalism. Nuclear fragmentation cross sections are calculated with a semi-empirical abrasion-ablation 
fragmentation model. The relation of the final computer code to the Boltzmann equation is discussed in the context 
of simplifying assumptions. 

The MCNPX code was developed at Los Alamos National Laboratory (LANL) and is described at their 
website 11 as “MCNP is a general-purpose Monte Carlo N-Particle code that can be used for neutron, photon, 
electron, or coupled neutron/photon/electron transport. Specific areas of application include, but are not limited to, 
radiation protection and dosimetry, radiation shielding, radiography, medical physics, nuclear criticality safety, 
detector design and analysis, nuclear oil well logging, accelerator target design, fission and fusion reactor design, 
decontamination and decommissioning. The code treats an arbitrary three-dimensional configuration of materials in 
geometric cells bounded by first- and second-degree surfaces and fourth-degree elliptical tori. Point-wise cross- 
section data typically are used, although group-wise data also are available. For neutrons, all reactions given in a 
particular cross-section evaluation (such as ENDF/B-VI) are accounted for. Thermal neutrons are described by both 
the free gas and S(alpha,beta) models. For photons, the code accounts for incoherent and coherent scattering, the 
possibility of fluorescent emission after photoelectric absorption, absorption in pair production with local emission 
of annihilation radiation, and bremsstrahlung. A continuous-slowing-down model is used for electron transport that 
includes positrons, x-rays, and bremsstrahlung but does not include external or self-induced fields. Important 
standard features that make MCNP very versatile and easy to use include a powerful general source, criticality 
source, and surface source; both geometry and output tally plotters; a rich collection of variance reduction 
techniques; a flexible tally structure; and an extensive collection of cross-section data. MCNP contains numerous 
flexible tallies: surface current & flux, volume flux (track length), point or ring detectors, particle heating, fission 
heating, pulse height tally for energy or charge deposition, mesh tallies, and radiography tallies.” 
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III. Results 


A. Two-layer Configurations 

There were three arrangements investigated for the two-layer configurations. The outer layer was either 
aluminum (Al), Tantalum (Ta), or Tungsten (W) and corresponded to 150 mils of thickness. The inner layers were 
high density polyethylene (HDPE) at 5 or 10 g/cm 2 thickness. These layers were compared with the same 
thicknesses of Al or HDPE. 

The results for the Al-HDPE system (Figure 2) show that the single layer aluminum case has the worst dose, 
while the Graded-Z Al-HDPE and single layer HDPE systems are comparable in dose. 
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Figure 2. Al-HDPE two-layer system compared with Aluminum and HDPE at the same thicknesses. 

The results of the Ta-HDPE system (Figure 3) show similar results to the Al-HDPE case. For the GCR 
environments, the Ta-HDPE system with 5 g/cm 2 of HDPE is slightly superior to the single layer of HDPE. 
However, the Ta-HDPE case with 10 g/cm 2 of HDPE is comparable to the single layer of HDPE. For the SPE 
environment, the reverse is true. The Ta-HDPE system is now comparable with the single layer of aluminum and 
the single layer of HDPE is superior for dose mitigation. 
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Figure 3. Ta-HDPE two-layer system compared with Aluminum and HDPE at the same thicknesses. 


The results of the W-HDPE system (Figure 4) show similar results to the Ta-HDPE case. For the GCR 
environments, the W-HDPE system with 5 g/cm 2 of HDPE is slightly superior, but the W-HDPE system with 10 
g/cm 2 of HDPE is more comparable to the single layer of HDPE. For the SPE environment, the W-HDPE system is 
more comparable to the single layer of aluminum, and the single layer of HDPE is the best radiation mitigator. 
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Figure 4. W-HDPE two-layer system compared with Aluminum and HDPE at the same thicknesses. 

The two layer configurations discussed in this section are similar to historical spacecraft thicknesses for radiation 
shielding purposes. The next section will consider very large thicknesses of water and whether the larger 
thicknesses change these results. 


B. Very Large Thickness Configurations 

For the very large thickness configurations, there were six designs analyzed that consisted of combinations of 
water with a borated polyethylene known as SWX. The three thicknesses of water were 100, 200, and 300 g/cm 2 
and the two thicknesses of SWX considered were 10 and 20 g/cm 2 . For comparison, a single layer of water at the 
same overall thicknesses was also considered. The results (Figure 5) show that the FfO-SWX systems are 
comparable to the single layer of H 2 0 at the same thicknesses. 
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Figure 5. H20-SWX very large thickness system compared with single layer of H20 at the same thicknesses. 


C. Verification and Validation 

For the verification and validation, high energy proton particles produced by SPE events and the proton 
component of the GCR are transported through shield material and tissue equivalent geometry using the Monte 
Carlo nuclear particle transport computer code, MCNPX. The MCNPX simulations were performed for the two- 
layer configurations and one case of the very large thickness configuration, for each environment case. 

The radiation source (GCR or SPE) spectra is modeled as an uniformly distributed planar source propagating 
from left to right along z axis toward a simple cylinder (of 15 cm radius). The input spectra source ranges from 0.01 
MeV to 2 GeV. The radial surface of the cylinder is modeled as a reflecting boundary to simulate an infinite material 
in the radial direction and to accomplish a one dimensional simulation. This allows for comparisons between 
MCNPX and HZETRN results. The MCNPX simulations are performed with the following particles transported 
explicitly: neutrons, protons, photons, electrons and light particles including deuteron, triton, helium-3, and alphas. 
The photon-electron transport is not present in the HZETRN version used in this study. 

The MCNPX code uses the mix-match algorithm that utilizes nuclear data tables to transport particles at certain 
energy ranges, and intranuclear cascade (INC) physics models to handle interactions at high energy regions. 12 The 
combinations of flux tally (F4), energy depositions per track length tally (F6), and cylindrical mesh tally were used 
to generate results. The F4 tallies were used to evaluate energy spectra at a transition point between shielding, at 
tissue phantom surfaces, and at various phantom depths. The F6 tally is used to calculate absorbed dose, averaged 
over a cell volume, resulting in dose (Gy/source particle/s). 

All calculations presented in this study are performed with MCNPX version 2.7.0. All simulations were 
performed at a total of 2 million histories. The calculated results were set to meet a minimum of 10% or less 
statistical uncertainty. The corresponding error bars were plotted in the curves, but in most cases, they were not easy 
to see because of its trivial values. 

Figure 6 shows the plot of neutron flux (secondary neutrons produced, resulting from the primary protons of the 
SPE) for the five different shielding cases. It is evident that the neutron fluxes (10-150 MeV) in the water shields are 
substantially lower as compared to the Graded-Z shield cases. Figure 7 and Figure 8 show the flux plots for the 
GCR (solar minimum and maximum) cases, respectively. Again notice the neutron fluxes for the water shield to be 
lower compared to the Graded-Z shield cases. Additionally, we notice (more pronounced) that the SWX layers in 
combination with water produce a higher thermal neutron flux as compared to the water case. This may be due to a 
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higher primary proton interaction with the shield material, resulting in a higher secondary neutron production. When 
these results are viewed in context with the dose results, it shows that the does is lower for the SWX cases when 
compared to just the water case. This suggests that the more damaging protons are removed by the SWX, providing 
additional shielding. 
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Figure 6. October 1989 SPE Band fit neutron flux. 
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Figure 7. 1977 GCR (solar minimum) neutron flux 
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Figure 8. 1982 GCR (solar maximum) neutron flux. 

Table 1 shows the summary of MCNPX dose results for the SPE. Clearly the dose values show a consistent 
trend to the HZETRN results. Overall, the MCNPX dose values are 15% to 30% higher than the HZETRN dose 
values. These results are consistent with previously published work. 12, 13 The absence of EM transport in HZETRN 
is one of the major reasons for these differences. Also the 2005 version of HZETRN does not include the upgrades 
to the neutron and light ion transport now available in the 2010 release. 14 This also contributes to the differences 
between the results. When the ratios of the relative thicknesses of the shields are compared, the difference between 
MCNPX and HZETRN is much smaller (5-10%). This is stated as the D ratio. This difference may be due to the 
differences in the transport calculation algorithms and in the cross-sections used by the two codes. The differences 
are higher for the 1982 GCR solar maximum results (Table 3), approximately 50% higher. Further investigation is 
necessary to understand why the difference is much greater than what we observe for the SPE and the 1977 GCR 
solar minimum results. 
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Table 1. Dose values for October 1989 SPE Band fit. 


1.0287 g/cm 2 Al 

HZETRN 

Dose, cGy D ratio 

MCNPX 

Dose, cGy D ratio 

Difference 

Dose D ratio** 




5 g/cm 2 PE 
10 g/cm 2 PE 

114.20 0.40 

45.94 

135.71 0.45 

60.44 

16% 10% 
24% 

6.3237 g/cm 2 Ta 




5 g/cm 2 PE 
10 g/cm 2 PE 

71.17 0.47 

33.41 

83.36 0.50 

41.83 

15% 6% 

20% 

7.3343 g/cm 2 W 




5 g/cm 2 PE 
10 g/cm 2 PE 

65.64 0.48 

31.51 

77.50 0.51 

39.42 

15% 6% 

20% 

100 g/cm 2 H20 




0 g/cm 2 SWX* 
10 g/cm 2 SWX* 
20 g/cm 2 SWX* 

0.74 

0.59 0.81 

0.49 0.67 

0.56 

0.42 0.75 

0.33 0.59 

24% 

29% 7% 

33% 12% 


*here the HZETRN dose is higher than MCNPX 

**the D ratio is the ratio of dose values from 10g/cm2 shield to the 5 g/cm2 


Table 2. Dose values for 1977 GCR (solar minimum). 


HZETRN 

Dose, cGy/day D ratio 

MCNPX 

Dose, cGy/hr D ratio 

Difference 

Dose D ratio 

1.0287 g/cm 2 Al 



5 g/cm 2 PE 0.0243 1.021 

10 g/cm 2 PE 0.0248 

0.03206 1.017 

0.03261 

24% 0% 

24% 

6.3237 g/cm 2 Ta 



5 g/cm 2 PE 0.0245 1.011 

10 g/cm 2 PE 0.0247 

0.03223 1.003 

0.03233 

24% -1% 

24% 

7.3343 g/cm 2 W 



5 g/cm 2 PE 00245 1009 

10 g/cm 2 PE 0.0247 

0.03227 1.000 

0.03228 

24% -1% 

24% 

100 g/cm 2 H20 



0 g/cm 2 SWX 0.0155 

10 g/cm 2 SWX 0.0145 0.932 

20 g/cm 2 SWX 0.0134 0.864 

0.01662 

0.01523 0.917 

0.01393 0.838 

7% 

5% -2% 

4% -3% 
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Table 3. Dose values for 1982 GCR (solar maximum). 


1.0287 g/cm 2 Al 

HZETRN 

Dose, cGy/day D ratio 

MCNPX 

Dose, cGy/day D ratio 

Difference 

Dose D ratio 




5 g/cm 2 PE 
10 g/cm 2 PE 

0.0031 1.0645 

0.0033 

0.0070 1.0570 

0,0074 

56% -1% 

56% 

6.3237 g/cm 2 Ta 




5 g/cm 2 PE 
10 g/cm 2 PE 

0.0032 1.0534 

0.0034 

0,0072 1.0373 

0,0075 

55% -2% 

55% 

7.3343 g/cm 2 W 




5 g/cm 2 PE 
10 g/cm 2 PE 

0.0032 1.0503 

0.0034 

0,0073 1.0330 

0.0075 

55% -2% 

55% 

100 g/cm 2 H20 




0 g/cm 2 SWX 
10 g/cm 2 SWX 
20 g/cm 2 SWX 

0.0029 

0.0028 0.9497 

0.0026 0.8960 

0.0053 

0.0050 0.9448 

0.0045 0.8651 

44% 

44% -1% 

42% -4% 


IV. Conclusions 

This study investigated the use of Graded-Z materials for radiation mitigation of proton environments. The 
results showed some slight dose reductions in using Tantalum or Tungsten with 5 g/cm 2 of HDPE in a GCR 
environment. However, the Graded-Z configuration was actually detrimental with an SPE environment when 
compared with a single layer of HDPE. In the very large thickness configurations, the results were comparable. 

The verification and validation showed similar trends between HZETRN and MCNPX. The differences in 
values between the codes were due to known discrepancies in methods of calculation. Further investigation will be 
needed to elucidate why there are larger differences for the GCR solar minimum environment case. 

Furthermore, since this study has only focused on the dose, additional work needs to be implemented to 
understand whether the incident radiation (protons) or the production of secondary radiation (neutrons) is 
responsible and how this will affect the dose equivalent. Furthermore, an analysis of the mass of these 
configurations is important because if we can acquire comparable radiation mitigation characteristics with a Graded- 
Z configuration, but decrease the overall mass of the system, there is still value in pursuing these Graded-Z 
structures. In addition, future work will also investigate three-layer systems to evaluate whether they provide 
advantages over two-layer or single-layer systems. 
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